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Process For Preparing Polyethylene 

5 BACKGROUND OF THE INVENTION 

1. Field of the Invention 

[0001] The present invention relates to a process for preparing polyethylene. 
The present invention particularly relates to preparing high density polyethylene 
using a slurry loop reactor. 

10 

2. Background of the Art 

[0002]lt is well known in the art of manufacturing chemicals in general and 
polymers in particular to use a reactor consisting of a pipe containing a flowing 
stream of reactants. One such pipe reactor is known as the slurry loop reactor. 
15 In a slurry loop reactor, reactants are feed into a pipe containing a solvent and a 
catalyst. The admixture of solvent, reactants, and catalyst are continuously 
recycled through the pipe, hence the term "loop," with the reaction product being 
continuously removed. 

[0003] In the manufacture of high density polyethylene, particles of polyethylene 
20 can be removed from the reactor by means of a settling leg. U. S. Patent No. 
3.293,000 to Marwil discloses the use of settling legs in a loop reactor for 
manufacturing an ethylene butane copolymer. In this reference, the settling leg 
functions to allow a portion of the polymer slurry to escape from the loop into 
outlets that allow the particles to be gravimetrically removed. 
25 [0004] While technology for preparing polymers has changed with respect to 
catalysts and reactants, the same general loop reactor technology employed in 
the 1 960s is still in use. For example, U.S. Patent No. 4,61 3,484 to Ayres, et al., 
discloses an improvement to settling legs, also referred to as accumulator legs, 
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is disclosed. This technology is still generally current and in present use, 
particularly for the production high density polyethylene. 
[0005] Despite loop reactors having been used for a long time, the process of 
starting a running a loop reactor is not without problems. U. S. Patent No. 
5 6,303,825 to Gerlich, et al. discloses a method of starting up a loop reactor 
system. This reference suggests using a microprocessor to automate starting 
up and running a loop reactor. One suggested algorithm recognizes the end of 
the startup of the reactor and then controls the reactor at steady state. 

10 SUMMARY OF THE INVENTION 

[0006] In one aspect, the present invention is a process for producing 
polyethylene using a slurry loop reactor comprising using a mathematical model 
to predict a plurality of process control parameters. The prediction is based on 
the desired product properties, reactor contents mass balance, and reactor 
15 characteristics. The slurry loop reactor is controlled using the predicted process 
control parameters. 

[0007] In another aspect, the present invention is a controller programmed with 
such a mathematical model for use in controlling a polyethylene slurry loop 
reactor. 

20 [0008] In still another aspect, the present invention is a process for optimizing 
the configuration of a polyethylene slurry loop reactor comprising using such a 
mathematical model to design the reactor. 

DETAILED DESCRIPTION OF INVENTION 

25 [0009] In one embodiment, the present invention is a process for producing 
polyethylene using a slurry loop reactor. The polyethylene produced can be 
high density polyethylene. The slurry loop reactors useful with the process of the 
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present invention include those using settling legs to remove polymer particles 
from the reactor. These reactors are also sometimes referred to as Phillips 
Reactors and manufacturing processes using such reactors are sometimes 
referred to as Phillips Processes. Loop reactors not having such legs can also 
5 be used with the method of the present invention. 

[0010] In the practice of the method of the present invention, a mathematical 
model is used. In one embodiment, this model uses mass balance of reactor 
contents, reactor geometry, catalyst kinetics, and settling phenomena in the 
outlet settling legs to determine relationships between inlet feed rates and 
1 0 reactor geometry on production conditions. The relationship between these can 
be defined using a series of equations having a plurality of variables. The 
variables that can be employed in a process of the present invention include: 





E, = 


ethylene flow into reactor 


15 


Eo = 


ethylene flow out of reactor 




// 


isobutane flow into reactor 




lo = 


isobutane flow out of reactor 




Hi = 


hexene flow into reactor 




H 0 = 


hexene flow out of reactor 


20 


H c = 


hexene conversion in reactor 




hi = 


hydrogen flow into reactor 




h 0 = 


hydrogen flow out of reactor 




h c = 


hydrogen conversion in reactor 




Po = 


polymer flow out of reactor 


25 


Lo = 


liquid flow out of reactor 




Tj w = 


total mass flow into reactor 




Tow — 


total mass flow out of reactor 




T 0 V = 


total volume flow out of reactor 




Ci 


catalyst flow into reactor 


30 


C 0 


catalyst flow out of reactor 




C E = 


concentration of ethylene in the reactor liquid 




C H = 


concentration of hexene in the reactor liquid 




C h = 


concentration of hydrogen in the reactor liquid 




T 


temperature of reactor 
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P 


pressure of reactor 




S w = 


weight concentration of solids in the reactor slurry 




• s v = 


volume concentration of solids in the reactor slurry 




S ws — 


weight concentration of solids in the settling leg 


5 




solids bed 




N, = 


number of settling legs 




V 


reactor volume 




d, = 


settling leg diameter 




h, = 


settling leg height 
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r b 


bulk density of reactor polymer 




r p 


density of reactor polymer 






density of reactor liquid 




r s 


density of reactor slurry 




t s 


residence time of reactor solids 


15 


k 


catalyst activity 




P 


catalyst productivity 




d c = 


catalyst diameter 




fa = 


catalyst feed factor 




k c = 


catalyst activity factor 


20 




terminal velocity of settling polymer 




Ps = 


polymer settling rate 




Pn = 


rate of polymer leaving reactor that is not part of 
the settling leg solids bed 




T s = 


rate of slurry leaving reactor that is part of the 


25 




settling leg solids bed 




T n = 


rate of slurry leaving reactor that is not part of the 
settling leg solids bed 




m, 


viscosity of reactor liquid 




r c 


density of catalyst 


30 


dp = 


polymer diameter 




N Ar = 


Archimedes number for polymer settling in settling 
leg 




N Re = 


Reynolds number for polymer settling in settling 
leg 


35 




acceleration due to gravity 






cross sectional area of a settling leg occupied by 
polymer 




A = 


cross sectional area of a settling leg 
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[0011] The mathematical model of the present invention is developed based 
upon these variables. For example, a given slurry loop reactor can have the 
following known values: 



p _ 


OOUOH IU/iI 


1 


COQfifi ih/h 
?C90u lU/l 1 


U£ = 


O.UOo/ WI /o 


Uh = 


u.o4oyo WT /o 


Uh - 


u.uu/oo iy wi/o 


U - 

n c - 


71 


H c = 


. 50% 


7 


210 F = 98.9 C 


P 


600 psig = 41 .8 atm 




37.668 wt% 


V 


100 m 3 


d, = 


10 in 


h, = 


15ft 


r b 


0.45 g/cc 


b 


0.89 g/cc 


P 


1624.9 lb HDPE/lb cat 


do = 


100 m 


r c 


2.2 g/cc 


0 


9.80 m/s 2 



25 

[0012] Relationships between the variables can be calculated as follows. 
Assuming mass balance around the reactor, 

T -(l +p Wo 

* tw \ o * o / v 

30 

T iw =L 0 +P o = 52966 lb/h + 58034 lb/h =11 1000 lb/h 

The component mass flows in the outlet liquid can be calculated based on 
known outlet concentrations. Hydrogen flow is negligible portion of total liquid 
35 flow. 
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E 0 



10 



100 


= — — — 

100 


L 0 C H 


(52966 lb/hX0.54695)_ 


100 


100 


L 0 C„ 


(52966 lb/hX0.0075519) 


100 


100 




I 0 +E 0 +H 0 =L 0 



= 41b/h 



I =L -E -H = 52966 lb/h- 2677 lb/h- 290 lb/h = 50000 lb/h 

O 0 o o 

[0013] Assuming mass balance of isobutane around the reactor, 



15 /, =I 0 =50000 lb/h 

[0014] Assuming constant conversion of hexene and hydrogen into polymer 
across the intended operating region, the hexene and hydrogen feeds into the 
reactor can be calculated: 

20 

= H 1 m-c M ) 

100 

IQOg, = (l00X2901b/h) = 100()lb/h 
' 100 -C H 100-71.03 

25 ° 100 

, = _iopA_ = (iooX4ib/h) =81b/h 

' 100 -C„ 100-50 
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[0015] Assuming that HDPE was produced from only the monomers ethylene 
and hexene, the ethylene feed into the reactor can be calculated from mass 
balance of ethylene around the reactor. 

E, = E 0 + [P a - (H i -H 0 )]= 2677 lb/h + [58034 lb/h - (1000 lb/h - 290 lb/h)] = 60000 lb/h 



10 [0016] Using the catalyst productivity and catalyst mass balance around the 
reactor, the catalyst feed to the reactor can be calculated. 



Co 



15 c 58034 lb/h =35.7151b/h 

p 1624.9 lb/lb 

c,-c 0 =0 

c . =c o =35.715 lb/h 



20 



[0017] An empirical correlation of isobutane liquid density, developed by Lewalle 
can be use with the process of the present invention. This correlation expresses 
liquid density (g/cc) as a function of liquid compositions (wt%), temperature (C), 
and pressure (atm) in the reactor: 



25 
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10 



p, = 0.5767 + 2.273e - 04P - C E 



3.991e-04 + 1.126e-04P 

- 2.496e - 06P 2 + 1.679e - 08P 3 



-rc t 



f4.943e-06 + 5.33e-06P 
-1.714e-07P 2 +1.307e-09P 3 
7.146e-05-5.734e-06P 
+ 1.254e - 07P 2 - 8.388e - 10P 2 



- P(7.718e - 04 + 2.362e - 06P) 

'8.449e-07-3.265e-07P ^ 



-TC f 



+ 9.876e - 09P 2 - 7.522e - 1 1P : 



- P 2 (7.646e - 06 - 6.735e - 08P)+ T 2 C t 

1.435e-08-4.099e-09P 
+ 1.197e-10P 2 -9.003e-13P 3 



r3.932e-07-5.549e-08P 
+ 1.269e-09P 2 -8.579e-12P 3 



+ T 2 C 2 



+ 0.001C, 



( 



0.001 + 



1- 



v 



P-15 
5C F +60 



( T \ 



2.5- 



[0018] The liquid density, upon substitution of variables, simplifies to 

p t =0.431 g/cc 

[001 9] Using the calculated liquid density, the assumed polymer density, and the 
solids concentration, the slurry density can be calculated: 



lOOp.p, 



(l00X0-89g/ccX0-431g/cc) 



Ps 100yO,+S wV 9 ( -pJ (l00X0.89g/cc)+(37.668X0.431g/cc-0.89g/cc) 

p s =0.535 g/cc 

15 [0020] The volumetric flow rate of slurry out of the reactor can be calculated 
from the mass flow rate and density of each component: 
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T — — 4- — 

ov ~~ n ■ ~ 

P P . Pi 



f 58034 lb/h | 5296611 
0.89 g/cc 0.43 1 g/cc J[ 





f Eal 1 


[ h 1 




^3785cc ; 


^60 min ; 



= 376 gpm 



5 [0021 ] The volumetric solids concentration in the reactor is: 



10 



s., =■ 



100 



100 



1+^ 

Pi 



100 



-1 



1 | ( 0.89 g/cc Y 100 l 
+ [ 0.431 g/cc y 37.67 > 



= 22.64 vol% 



[0022] The solid residence time is calculated as follows: 

Vp s S w _ (l00m 3 )(0.535g/ccX37.668) le06cc Jb 60_min _ 4< . Q 



100/> 



(100X58034 lb/h) 



m 3 453.6 g h 



nun 



[0023] Assuming HDPE production kinetics as being first-order for ethylene 
15 concentration in the liquid and catalyst flow rate, the kinetic rate can be 
expressed as: 



P„=kc 0 C E t s 



20 [0024] Substituting the catalyst productivity and rearranging the kinetic rate 
equation, the catalyst activity is calculated as follows: 



Co 



25 



k = 



1624.9 lb/lb 



(45.9minX5.0537 wt%) 



= 7.00 lb/lb/min/wt% 
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[0025] The catalyst feed factor is a function of other measured reactor variables 
that makes it independent of catalyst type or catalyst activity. It is expressed as: 



5 



10 



15 



[0026] Using mass balance of catalyst around the reactor and rearranging the 
kinetic rate equation yields the following equation for catalyst feed factor: 

r , , P o 580341b/h . , m 

f = kc, = kc„ = — — = 7 w r = 250 lb/h/min/wt% 

0 t s C e (45.9minX5.0537wt%) 

[0027] The catalyst activity factor quantifies the reactivity of the specific catalyst 
with respect to its size. It is expressed as: 



k 7.001b/lb/min/wt% le061b _„ . , M , 3 

k= — r = = 7.00 lb/MMlb/min/wt%/ i u 

c (lOO^) 3 IMMlb 



[0028] The slurry leaving the settling leg consists of two portions: the settled bed 
of polymer and liquid in the bottom of the settling leg, which is a larger 
concentration of solids than the reactor, and the polymer and liquid above this 
20 bed, which is the same concentration of solids as in the reactor. The 
concentration of polymer in the settling bed solids bed is a function of the 
polymer and liquid densities as well as the polymer bulk density and is 
calculated as follows: 



100 100 wt% „ 
25 S ws = ^ = r = 67.9 wt% 



1 + P, 



1 1 

[P b ~P P ) 



1 + (0.431 g/cc] 



0.45 g/cc 0.89 g/cc 
[0029] A mass balance of the two portions of the outlet slurry is 
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T = T +T 

[0030] A mass balance of the two portions of polymer leaving in the outlet slurry 

5 is: 

p - p + p 

1 O S ft 

[0031] Using the known polymer concentrations of these two portions and 
10 solving the set of simultaneous equations, the amount of polymer leaving the 
settling leg from the settled bed of polymer can be calculated. This flow rate also 
indicates the polymer settling rate in the settling leg: 

100P 

T. = 



15 



s S 

ws 

100P. 



T. =■ 



100P loop 



T = - + ■ 

x ow 



n 



^ 100P, lOOfo-Pj 
20 T ow =—*- + 



aw 

' ' ws 



T ow P 0 1110001b/h 58034 lb/h 

100 S w = 37.7 = 3646Q lb/h 

5 J l_ _1 1_ 

S ws S w 67.9 37.7 

[0032] The viscosity of isobutane liquid as a function of temperature is: 



25 
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,u,(cP) = e 



-7.3891+- 



2582.6 
' 172.23+T(K) _ e 



-7.3891+ 



2582.6 



172.23+98.9+273.15 = 0.0711cP 



[0033] An equation offered by Mignon to describe the relation between the size 
of the growing polymer particle, the catalyst particle size, and the catalyst 
productivity can be used with present invention. The polymer diameter is 
expressed as: 



10 



15 



d p =0A2d c 



( V' 3 

PcP 



= (0.42X100// 



(2.2g/ccXl625)" 



0.89 g/cc 



1/3 



= 666// 



[0034] Assuming the polymer settling in the settling leg to be steady-state and 
the polymer particles are close to spheres, the settling rate can be described by 
the following equations. The Archimedes number for the polymer settling 
through the isobutane liquid is: 



N Ar = 



d 3 p lg (p p -p,) 



= (666 /if (0.43 1 g/cc)(9.80 m/s 2 g/cc - 0.43 1 g/cc) ( cP 



(0.071 lcP) 2 



V Mm ^ 



g/m/s 



le06^ 



'100 cm ^ 6 



lm 



20 



N M =1.13e05 



[0035] An equation by Dallavalle expresses the drag coefficient for this type of 
flow up to Reynolds number of 2e 05 . This equation can be used to calculate the 
Reynolds number from the Archimedes number: 



12 
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iV Re = ^14.42 + 1.827^ -3.798 J = + 1.827>/l.l3e05 -3.798 



AT Re =453 



5 [0036] The Reynolds number confirms that the above equation is valid for this 
flow regime. The velocity of the settling polymer can now be calculated using the 
Reynolds number: 



N -^L 

iV Re - 



Ml 

10 



v, = 



N Re /i, (453X0-07 llcP) g/m/s le06/< 



m 



100 cm 



3.28 ft 



m 



= 0.369 ft/s 



' d p p l (666//X0- 43 lg/cc) cP m 

[0037] The expression relating the mass flow of settling solids to the velocity of 
those solids is: 

15 - : 

P s =N,A lp v t p p 

[0038] The area of the settling leg occupied by polymer is: 



20 A <> = TbF = ^' X 

[0039] Substituting this area into the mass flow of settling solids equation yields: 



p * = ik N ' d ' Xv ' p > 

25 

[0040] Rearranging this equation, the number of settling legs required for this 



reactor geometry and these production conditions is: 
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7id?S v v t p p 

(400X36460 lb/h) 453.6g h ft ( in V = 4 

, ~^-(l0in) 2 (22.6)(0.369ft/sX0.89g/cc)' lb '3600s 12 in [ 2.54cm J 

[0041] In the practice of the process of the present invention, these equations 
5 can be used to solve for numerous combinations of knowns and unknowns. For 
example, if the number of settling legs is known, then the diameter of the settling 
legs can be calculated. This is a simple calculation using the same sequence of 
equations as above. 

[0042] For other variables, the equations above may be used in alternative 
l o sequences. An unknown variable, for example the ethylene concentration in the 
liquid, must first have an estimated value and then the sequence used in an 
iterative fashion until the ethylene concentration converges on a solution that 
satisfies all the equations. This type of calculation can be performed using the 
equations detailed above with an automated spreadsheet. 
15 [0043] In one embodiment, the present invention is a controller for an industrial 
high density polyethylene slurry loop reactor. While the present invention can be 
used in manual or spreadsheet form, it can also be incorporated into the logic 
circuits of a controller. In one embodiment, the controller is a neural net or other 
artificial intelligence (Al) controller. When using the model of the present 
20 invention with a spreadsheet, it can be desirable to incorporate additional 
functionality into the spreadsheet using, for example, VISUAL BASIC®, or some 
other compatible computer program. It is also within the scope of the present 
invention that the model be incorporated into a computer program such as 
Fortran or C++. 

25 
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[0044] Many operations in a chemical process are routinely controlled using a 
Proportional. Integral Derivitive (PID) controller. These controllers, while very 
useful in a conventional setting, are not preferred for the present invention. 
Such controllers are designed and used to monitor a single process variable and 
5 perform their control function based upon that single variable. In an 
embodiment of the present invention, many PID controllers, are used in 
conjunction with a second controller that can receive data from the PID controller 
and then reprogram the PID controller based upon the total inputs from the 
reactor sensors and controllers. 

10 [0045] In another embodiment, the present invention is implement using an Al 
controller capable of accepting multiple inputs and sending multiple outputs. 
Preferably, the controller is a neural network model based controller. For 
example, one such controller is a controller using Process Perfecter® software 
developed by Pavilion Technologies. 

15 [0046] While the automated controllers can be desirable for their quick . 
response, in effect controlling the reactor in .real time, the use of the present 
invention manually should not be discounted. When incorporated into a 
spreadsheet, the process of the present invention for controlling a reactor can be 
very useful, particularly when changing polymer grades or production rates. In 

20 either case, the process of the present invention can be used to optimize and 
control a loop reactor to produce polyethylene with desirable properties. In an 
alternative embodiment, the model can also be used to design a loop reactor 
such that settler size, cycle time, reactor volume and energy consumption can all 
be optimized. 

2 5 [0047] While the primary subject of this application is directed to the production 
of polyethylene, it is well known that loop reactors can be used to produce other 
types of polymers. For example, polypropylene can also be produced using the 

15 
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method of the present invention and is within the scope of the present invention. 
Any polymer that can be produced using a loop reactor can be produced using 
the method of the present invention. 

5 EXAMPLES 

[0048] The following examples are provided to illustrate the present invention. 
The examples are not intended to limit the scope of the present invention and 
should not be so interpreted. Amounts are in weight parts or weight 
percentages unless otherwise indicated. 

10 

EXAMPLE 1 

[0049] The mathematical model described above is placed into a spreadsheet 
and automated to allow for quick calculation of solutions. The intended 
operating range of the reactor is used to form the ranges of the known variables. 

15 Within these ranges, the variables are randomized and the spreadsheet used to 
calculate solutions to several different combinations of production conditions. 
These results are placed into a large dataset and fed into the control software 
for the building of a neural network control model. The resulting neural network 
control model then possesses all of the inherent relationships between variables 

2 0 of the original equation-based model. 

[0050] The accuracy of the equation-based model is proven when the control 
model is placed online for closed loop control of key reactor compositions and 
rates ethylene concentration, percent solids, hexene concentration, hydrogen 
concentration, hexene/ethylene concentration ratio, hydrogen/ethylene 

25 concentration ratio and production rate. With the model-built controller, the 
variability of these control variables is decreased by 70% in comparison to 
conventional PID controllers replaced by the controller of the present invention. 
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[0051] The robustness of the model is demonstrated by the ability of the 
controller to perform successfully at from about full reactor capacity down to 
about half capacity and across a large range of compositions. The generic 
5 nature of the model with respect to catalyst using catalyst feed factor and 
catalyst activity factor to describe catalyst instead of catalyst feed flow, catalyst 
activity, and catalyst size is been shown by its ability to be used for Chromium, 
Ziegler-Natta, and Metallocene catalysts successfully with the same controller. 

10 EXAMPLE 2 

[0052] A loop reactor is modeled according to the method shown above, with 
relationships established between the variables. Assumptions are made and 
then confirmed upon experimental observation. A computer spreadsheet is 
programmed with the model. The interface for this model is displayed in 
15 APPENDIX A as the RX Model. The formulas are displayed in APPENDIX B. 

EXAMPLE 3 

[0053] A computer spread sheet is programmed as in Example 2 except that it is 
set up to allow for the input of feeds to the reactor to calculate the outlet flow 
2 0 rates and compositions. The interface for this model is displayed in the 
APPENDIX A as the RX Model Iterative. The formulas are displayed in 
APPENDIX C. 

25 
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